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Abstract 

The amplitudes for the pp — > ppir process applicable at high energy are discussed in detail. 
Both diffractive bremsstrahlung (Drell-Hiida-Deck-type model), photon-photon and photon-omega 
exchange mechanisms are included in the calculation. Large cross sections of the order of mb 
are predicted. The corresponding differential cross sections in rapidities and transverse momenta 
of outgoing protons and pions as well as relative azimuthal angle between outgoing protons are 
calculated for RHIC and LHC energies. The hadronic bremsstrahlung contributions dominate at 
large (forward, backward) pion rapidities. The diffractive nonresonant background contributes 
at small tt°p invariant mass and could be therefore misinterpreted as the Roper resonance. We 
predict strong dependence of the slope in t (squared four-momentum transfer between ingoing 
and outgoing proton) on the mass of the supplementary excited 7r°p system. At high energies 
and midrapidities, the photon-photon contribution dominates over the diffractive components, 
however, the corresponding cross section is rather small. The photon-odderon and odderon-photon 
contributions are included in addition and first estimates (upper limits) of their contributions to 
rapidity and transverse momentum distribution of neutral pions are presented. We suggest a search 
for the odderon contribution at midrapidity and at p± n o ~ 0.5 GeV. Our predictions are ready 
for verification at RHIC and LHC. The bremsstrahlung mechanisms discussed here contribute also 
to the pp — > p(mr + ) reaction. Both channels give a sizable contribution to the low-mass single 
diffractive cross section and must be included in extrapolating the measured experimental single 
diffractive cross section. 
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I. INTRODUCTION 



In the past we have studied exclusive production of pairs of charged pions in the pp — » 
ppiT + 7r~ [l|, [2| and pp — > nnir + ir + ^ processes. There pomeron-pomeron or pomeron- 
reggeon exchanges are the dominant mechanisms. The corresponding cross sections for high 
energies are rather large, of the order of 100 fib. Such processes could be measured with 
the help of the main ATLAS or CMS detectors (for charged pions detection) , ALFA |ij or 
TOTEM Q detectors (for protons tagging), and the zero degree calorimeters (ZDCs) (for 
neutrons detection). Furthermore, the proposed forward shower counters, to detect and 
trigger on rapidity gaps in diffractive events, would improve the measurements at the LHC 
significantly [6j. 

Here we consider another simple final state with one pion only. The exclusive process 
pp — > pp7r° was measured in detail only near to the pion threshold at the IUCF (Bloom- 
ington) Q, CELSIUS (Uppsala) @ and the COSY (Jiilich) [9]. The total cross section for 
single pion production grows from threshold to about 10 /A> at the cm. energy 3 GeV. 

Although only a few partial waves are involved close to threshold, the theoretical description 
is not easy (see e.g. [lOj and references therein). For a summary on near-to-threshold meson 
production experiments see (llj . 

What happens when the energy increases? In the range of center-of-mass energies \/s = 
3—10 GeV the nucleon resonances can be excited via meson exchange processes. Evidence of 
proton excitation can be observed in the pir mass spectrum (A + or N* + ). A nice summary 
of the intermediate energy data for pp — > ppir can be found in Ref.pjj]. In this region of 
energy the corresponding cross section systematically decreases which is consistent with the 
meson exchange picture. When energy increases further the role of many of the nucleon 
resonances diminishes and the mechanism becomes simpler. 

In Refs [l3j, llj] a study of pseudoscalar mesons produced centrally by the CERN-WA102 
Collaboration at ^/i = 29.1 GeV was performed. The results show that the 77 and 77' 
mesons appear to have a similar production mechanism which considerably differs from 
that for the n° production [l3[. To our knowledge this was never explained theoretically. 
The WA102 Collaboration concentrated on very central production of mesons and therefore 
measured protons with large Feynman xf = 2p\\/ yfs. This condition eliminates contribution 
of the diffractive mechanisms discussed in our paper. Reactions of this type pp — > pMp are 
expected to be mediated by double exchange processes, with a mixture of pomeron-pomeron, 
reggeon-pomeron, and reggeon-reggeon exchange. For instance, the 77 and 7/ mesons are 



produced dominantly by double pomeron exchange (see 15[ and references therein). For 
the central exclusive 7r° production at intermediate energies the p-u exchange may be the 
dominant mechanism. The p-a.2 exchange could be another potential candidate. The validity 
of these exchanges could be justified experimentally by the COMPASS Collaboration (see 
e.g., Ref.Q). 

In the present paper we wish to concentrate on the production of single neutral pions in 
the pp — > ppn° reaction at large energies (RHIC, LHC). We hope that this process could 
be measured, at least in some corners of the phase space, at the LHC. We shall refer also 
to pp — > p(mr + ) and np — > (pft~)p reactions measured at lower energies at Intersecting 
Storage Rings (ISR) and Fermilab in the 1970's (l7l-[2l| (for a nice review we refer to 22|). 
The mechanism of these reactions is closely related to the pp — > ppir reaction discussed 
here and will be therefore a good reference point for our calculation. As discussed in the 
past, the dominant hadronic bremsstrahlung-type mechanism is the Drell-Hiida-Deck (DHD) 
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mechanism for diffractive production of ttN final states in NN collisions 23] ; for a review, 
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see e.g. 22, 24]. The exclusive pion production mechanism is similar to pp — > ppuj [25| and 
pp —> ppj 



processes. 

The 7T° can be also produced by 77, 70;, and 7C exchanges. The soft odderon (O) couples 
very weakly to the nucleon. In Refs. [27], [28| the authors discussed some results of exclusive 
pseudoscalar meson production in high energy ep scattering. It was shown in [28[ (see also 
[29( 1 ) that odderon exchange leads to a much larger inelastic than elastic 7r° production cross 
section. As shown in Ref. [28], the photon exchange is larger than the odderon exchange 
only at very small transverse momenta of 7r°. In this paper we shall consider the odderon 
contribution in proton-proton collisions using a simple phenomenological approach for the 
odderon exchange. We shall discuss how it can be separated from the contribution of photon- 
photon fusion. 



II. THE AMPLITUDES FOR THE pp ppvr REACTION 



A. Diffractive bremsstrahlung mechanisms 
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FIG. 1. Diagrams of the 7r°-bremsstrahlung amplitudes driven by the pomeron exchange in proton- 
proton collisions: (a) pion exchange, (b) proton exchange, and (c) direct production. The direct- 
channel p* in (c) is an off-shell proton, not a proton resonance. Some kinematical variables are 
shown in addition. 



The diffractive bremsstrahlung mechanisms for the exclusive pp — > ppn° reaction are 
driven by pomeron (IP) as depicted in FigJT] and reggeon (M) exchanges. At high cm. en- 
ergies \/s the dominant contribution comes from the pomeron exchange. There are two 
processes when the 7r° meson emitted by one of the protons interacts with the second proton 
(diffractive tt° rescattering) , as depicted in Fig{T^a), and four processes in which protons 
interact and the n° emission may occur, see Figs{T](b) and[T](c). In general, the amplitudes 
of these processes may interfere but in the present case the interference is negligible as the 
two processes are well separated in rapidity as will be discussed in Section IHIl The Born 



1 Discussed here diffractive mechanisms of exclusive tt° production are similar to the diffractive mechanism 
of pp —> ppuj [25! and pp — ¥ ppj [2^] processes. 
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amplitudes for these processes, see Fig{TJ can be written as 

M X~X^fxT£° = U ( Pl > ^l)h5u(p a , Xa) S^fa) g^NN F^NN (ti)Fj^. w (ti) 

x (A^ N (s 23 ,t 2 ) + A^(s 23 ,t 2 )) /(2s 23 ) 

x (<?i + Ps)n «(P2, A 2 )7 m m(P6, Aft) , (2.1) 

M X n 7^t\!\onO = 9nNN «(pi, Xl)^ S N {jpl^i^u{jp a , K) F n NN*{Pli) F 1PN*n{Pu) 



x (A™( Sl2 ,t 2 )+A™(s 12 ,t 2 )) /(2s 12 ) 



x w(P2, A 2 )7 A1 u(pft, Aft) , (2.2) 
^£^A™t? = 9ttNNu(PI, Ai)i75<SW(Pif)7MPa, A a ) F nN * N (p\ f ) Fjp NN *(p\ f ) 
x (^( S ,i 2 )+^( S ,t 2 )) /(2a) 

x w(p 2 , A 2 )7 M M(pft, Aft) , (2.3) 

where A), u(p',X') = (p' , X')j° are the Dirac spinors of the incident and outgoing 
protons with the four-momentum p and the helicities A; normalized as u(p')u(p) = 2m p . The 
factors or ^ appear here as a consequence of using spinors. The four-momenta squared 

of intermediate particles are defined in FigHJand p\^ 2i = (p a ,b — P3) 2 , Pi/, 2 / = (pi,2 +P3) 2 , 
Q12 = (Pa,b — Pi,2) 2 , the four-momentum transfers along the pomeron line ti >2 = q\ 2 and 
s jj = (ft +Pj) 2 are squared invariant masses of the system. The propagators of the 
off-shell pion and proton are 

SAt) = 7 — " — 2 > $n(p 2 ) = l St^4 , (2.4) 



t — m 2 p 2 — mi 



where ^ = p M 7 M . 

The energy dependence of the elastic scattering A(s, t) was parametrized in the Regge-like 
form with pomeron {i = IP) and reggeon (i = M = f 2 , p, a 2 , u) exchanges, 

Af(s,t) = Vl Qs \ -J exp l-L-j , (2.5) 

where we use the scale parameter sq = 1 GeV 2 . In writing the above amplitudes (12. 5p we 
have omitted indices related to helicities as we have assumed helicity conservation in the 
rescattering process. If the energy in the np or the pp system is small, then the secondary 
trajectories are also important, e.g. we have in (12. ip term A 71 p = Ap + Af 2 . The strength 
parameters Cj, the values of signature factors % and the (linear) Regge trajectories «i(t) = 
«i(0) + a[t are taken from the Donnachie-Landshoff analysis [30j| of the total NN and np 
cross sections and are listed, e.g., in Ref.[l[. From the optical theorem we have a tot (s) = 
l/s!mA el (s,t = 0). The running slope for elastic scattering can be written as B(s) = 
Bf + 2a;in(s/soL where B% N = 9 GeV~ 2 , = 5.5 GeV" 2 and B% N = 6 GeV" 2 , 

B-g = 4 GeV" 2 [l|, |25| for pomeron and reggeon exchange, respectively. 



Usually a high-energy approximate formula is used in the literature in calculating differ- 



ential cross sections. We use a precise calculation of the phase space (see e.g., |3l|). This is 
important if one wants to go to lower energies and/or to large rapidities. As will be discussed 
in the next section, for this particular reaction the cross section has maximum just at large 
rapidities, where the often used formula is too approximate. In the high-energy limit we 
obtain 

(?i + Ps) M «(P2, X 2 )^u( Pb , Aft) ( qi + p 3 )^ (p 2 + p b Y5x 2 x b = 2s 23 6 X2 x b . (2.6) 
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In the bremsstrahlung processes discussed here the intermediate nucleons are off-mass 
shell. In the above equations the off-shell effects related to the non-point-like protons in the 
intermediate state are included by the following form factors: 

*V)= ,., ■ (2-7) 



(p 2 — m 2 ) 2 + A 



N 



Such a form was used e.g. in Ref.[32] for u photoproduction. In general, the cutoff parame- 
ters in the form factors are not known but could be fitted in the future to the (normalized) 
experimental data. From our general experience in hadronic physics we expect A at ~ 1 GeV. 
We shall discuss how the uncertainties of the form factors influence our final results (see Sec- 
tion nni). 



The pion-nucleon coupling constant g n NN is relatively well known [33j, l34|. In our cal- 
culations we take g\ NN jA. r K = 13.5. F 7T * NN (t) is a vertex form factor due to the extended 
nature of particles involved. Unfortunately, the off-shell form factor is not well known as it 
is due to nonperturbative effects related to the internal structure of the respective objects. 
This discussion of form factors applies also to the other JP7T*ir vertices. We parametrize 
these form factors in the simple exponential form, 

ft — m 2 \ 

F^ NN {t) = Fjp n , n (t) = exp f ) , (2.8) 

which is conventionally normalized to unity on the pion-mass shell and A w = 1 GeV is a 
reasonable choice. 

The pion exchange as a meson exchange is a correct description at rather low np energies 
while a reggezation of pion is required at higher energies. We propose to use a generalized 
pion propagator (see e.g. 35|) at an appropriate subsystem energy and t, 



S w (t) f3 M (s)S n (t) + p R ( 8 )V*(t, s) , (2.9) 
where the pion Regge propagator with the Euler's gamma function, 

[hS) 2rK(t) + l) sm(na n (t)) {s J ' ^ 

gives a suppression for large values of t. We have introduced extra phenomenological func- 
tions (3m (s) and (3r(s) to interpolate between meson and reggeon exchanges. We parametrize 
them as 

P M (s) = exp (-(s - s thr )/A 2 int ) , /3 fl (s) = 1 - (3 M {s) . (2.11) 

From our general experience in hadronic physics the parameter Aj nt = 1 GeV. The pion 
trajectory is taken as a n (t) = ot' n (t — m 2 ) with a slope parameter a' n = 0.7 GeV -2 . 

We improve the parametrization of p-exchange amplitude (12. 2p to reproduce the high- 
energy Regge dependence by the factor (si 3 /s t / ! , r )° !JV< ' Pll ' )_5 or by the factor (S23/ s thr ) aN( - P2i ^ 5 ( 
where the threshold factor s t hr = ( m P + m n o) 2 and the nucleon trajectory is «Ar(Pij 2 i) = 
-0.3 + a' N p 2 u 2i with a' N = 0.9 GeV -2 . 
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FIG. 2. Typical absorptive correction diagrams to: (a) pion exchange and (b) proton exchange. 
The stars attached to protons and ir° meson denote the fact that they are off-mass shell. 



B. Absorptive corrections 



Let us estimate absorptive corrections 5M. shown in Figj2j Diagrams which involve 
the elastic scattering of the incident protons are termed "initial-state" absorption. From 
physical reasons discussed in 22|, |36| the diagrams, when the transition of excited proton 
p* p7T° occurs inside, do not contribute significantly at high energies. Diagrams with 
the "final-state" absorption corrections provide the dominant absorptive effect [361] - In the 
quasieikonal approximation which takes into account contribution of elastic rescatterings 
the absorbed amplitudes can be expressed as 

AW-Pll, -P2±) = M(-p 1± , -p 2± ) - 5M(-p 1± , -p 2± ) , (2.12) 
where 5Ai for the diagrams with "initial-state" absorption is the sum of convolution integral 



r A ./initial state abs , 



-Pl±>~P2±, 



]ir 2 s 



d 2 k 



^a 



X 



M 
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A' a A£->AiA 2 7r°V Pl-L' P2±) 

and in the case of diagrams with "final-state" 

$Mi:f:Xx h M-pi±>-P2±) = 



M 
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\NN 
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(2.13) 



>AiA^ 



+ 



d 2 k 



Sl2 

_!_ KA (p — excnan S e ) | 
523 ^-^1^° 



~Pl±, -P2±) A K y 2 ^ XlX2 (s l2 , k ± ) 
~Pl±, -P2±) A lLx 2 (s23, k±) , 



(2.14) 



where the two-dimensional transverse vectors —pi± = ~Pi± + k± and — p 2± = —p 2 ± — k± are 
the transverse components of the momenta of final state protons and k± is the momentum 
transfer. A el (s,k±) is an elastic scattering amplitude given by Eq. (j2.5j) at an appropriate 
energy and for the momentum transfer k±. Since in our calculations we include effec- 
tive pomeron and reggeon exchanges, i.e. pomerons and reggeons describing approximately 
nucleon-nucleon or pion-nucleon elastic scattering, no explicit absorption corrections have 
to be included in addition. 

Experience from hadronic phenomenology (for several analyses of two-body reactions 



sec 



371 ]) suggest that the purely elastic rescattering taken into account by Eq. f)2.12p are 



insufficient, and inelastic intermediate states (screening corrections) lead to an enhancement 
of absorptive corrections. This is sometimes included in a phenomenological way by a factor 
A sc (X sc > 1). Taking into account absorption corrections, the DHD mechanism was shown to 
give a reasonable explanation for the main properties of the low-mass diffractive dissociation 
38j . The effect of the absorption in diffractive dissociation is also discussed in 39 . 
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C. 77 and "fu) exchange mechanisms 






FIG. 3. A sketch of the photon-photon (a) and photon-omega meson (b) exchanges induced 
production of 7r° in proton-proton collisions. 



In the following we wish to investigate competitive mechanisms to the diffractive mecha- 
nisms discussed in the previous subsection. The new mechanisms, never discussed so far in 
the literature, are shown schematically in Figj3j In the most general case the corresponding 
Born amplitudes read 



\A 77- exchange 
" /vl A a A 6 ^AiA 2 7rO 



eu{p u Ai)7 At w(p a , A a ) Fx(ti) 



X 



5V 



e 2 e^' pa qi , p q 2 ^F rr ^Ati,t2) 



x eu(p 2 , \ 2 )j v u(p b , X b ) Fi(t 2 ) , 



(2.15) 



x 



9w' 



fi V per 



(h,t 2 ) 



li t 2 — m u 

X g^NN U(p 2 , X 2 )l l 'u(p b , \ b ) F uNN (t 2 ) F(s 23 , t 2 ) , 

MV^l^l = g^NN u(pi, \i)^u(p a , A a ) F wNN {t!) J*(si 3 , ti) 



(2.16) 



X 



"5V 



x eu(p 2 , \ 2 )Yu(pb, At) Fi(i 2 ) 



(2.17) 



where the ~/*NN vertices are parametrized by the proton's Dirac electromagnetic form factor 



4m 2 p -2.79 1 



(4m 2 -t)(l-t/m 2 D ) 2 ' 



(2.18) 



where m 2 D = 0.71 GeV 2 is a phenomenological parameter. 

The central vertices involve off-shell particles. The t dependences of F 7 . 7 .^ 7r o(t 1 , t 2 ) 
electromagnetic off-shell form factor are the least known ingredients in formula ( I2.15p . It 
is known experimentally only for one virtual photon 77* — > 7r° (e.g., 40j). In the present 
calculation we use a vector meson dominance model inspired parametrization of the 7*7* — > 
7T° transition form factor, 



F^^^oiti, t 2 ) 



F rr ,o(0,0) 



(l-tt/mDil-h/ml) 



(2.19) 



where m p is the p meson mass. The form factor is normalized to F 7 * 7 * 7r o(0, 0) = 12 ^f > 
where N c = 3 is the number of quark colors and f n = 93 MeV is the pion decay constant. 

The coupling of the omega meson to the nucleon is described by the coupling constant 
9uNN /^tt = 10 an d the corresponding form factor is taken in the exponential form: 



9 

m 



lv uiNN 



F uNN (t) = exp ( -2-^ ) , (2.20) 



where A u nn = 1 GeV. The g^-y — 0.7 GeV 1 constant was obtained from the omega 



partial decay width as discussed in Ref. [25J. The •yu and wy form factors are taken in the 
following factorized form: 

™l (U- m 2 



F ruJ ^AtiM) = ~W ex P ^— -Si . (2.21) 



In practical calculations we take A W7r7 = 0.8 GeV [25[ as found from the fit to the yp — > up 
experimental data. 

At larger subsystem energies, s y - 3> s^ r , one should rather use reggeons than mesons. 
The "reggezation" is included here only approximately by a factor assuring asymptotically 
correct high energy dependence, 

/ <! \ I arCtan (( S ~ S ""')/ A thr)( Q JRW- 1 ) 

F(s,t)= — , (2.22) 

V Sthr J 

where A thr ~ 1 GeV, a^O) = 0.5, and a' R = 0.9 GeV~ 2 . 

In the high energy limit we can write a relatively simple formula of two-photon fusion 
amplitude squared and averaged over initial and summed over final spin polarizations (see 

0): 

I^^T = 4sV ^ F [^ 2) \F r ^Xh,t 2 )\ 2 | qi J 2 |q 2 J 2 sin 2 (0 12 ) , (2.23) 

where 4>\2 — 4>i ~ 02 is the azimuthal angle between the two outgoing protons. 

The amplitude for processes shown in Figj3] are calculated numerically for each point 
in the phase space. In calculating cross section we perform integration in log 10 (pu_) and 
log 10 (p2±) instead in pi± and p2±, which is useful numerically because of photon propagators. 



D. jO and O-f exchanges 

As will be shown in Section 1111} at the ir° midrapidity only the 77 — > 7r°, out of the 
mechanisms considered so far, contributes, i.e. the corresponding cross section is rather 
small. This gives a chance to search for yO and Oy exchange processes shown in FigJH 
The 7£> — y Tc°p reaction was proposed some time ago as a good candidate for identifying 
the odderon exchange, the C = — 1 partner of the pomeron They have predicted cross 



section of about 341 nb at the HERA energy. However, the search performed at HERA [45( 
was negative and found only an upper limit for this process o yp ^a v < 49 nb. Ewerz and 
Nachtmann [47] found an explanation of this discrepancy within a nonperturbative approach 
using approximate chiral symmetry and partially conserved axial vector current (PCAC). 
They have found that the amplitude for diffractive neutral pion production is proportional 
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to m 2 and vanishes in the chiral limit (m^ — y 0). They have estimated that the cross section 
is probably damped by a factor of 50 (see |46| ) compared to the early estimate in [28]. 
The exclusive production of neutral pions at midrapidity can be used to search for odderon 
exchange as well as to test the predictions of Ref . [47( . 




7r°(p 3 ) 




7r fe) 



FIG. 4. Diagrams with the photon-odderon (a) and odder on-photon (b) exchanges in the pp —y ppir 
reaction. 



The cross section for photon-odderon and odderon-photon exchanges can be estimated 
in the equivalent photon approximation similar as was done in Ref. 26( for other photon 
induced processes. In such an approach the distribution of the neutral pions can be written 

as 



-pi) 

-vl) , 



(2.24) 



where f(z) is an elastic photon flux in the proton; an explicit formula can be found e.g. in 
. In the formula above, zi/ 2 = ^=-exp(±?/) with m± = a/ + p\ . 
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The differential cross section jp — y 7r°p will be parametrized in the present paper as: 



da. 



7p— »7r u p 

~dt 



B 2 (-t)exp(Bt)(x. 



7p— >ir°p • 



(2.25) 



The differential cross section vanishes at t — which is due to helicity flip in the 7 — y 7r° 
transition. The slope parameter can be expected to be typically as for other soft processes 
B ~ 4 — 8 GeV~ 2 . At the LHC and at midrapidities typical energies in the photon-proton 
subsystems are similar as at the HERA. In the following we shall consider two scenarios: 
HERA upper limit (a lp ^. n o p = 49 nb) and Ewerz-Nachtmann estimate (a lp ^ n o p = 6 nb). 



III. RESULTS 

Now we present our calculations of cross sections and distributions of the exclusive ir° 
meson production in proton-proton collisions. The rapidity distributions of 7r° are shown in 
FigElat center-of-mass energies y/s = 45 GeV (ISR), 500 GeV (RHIC) and 14 TeV (LHC) 
for all processes considered in the present paper. We present results for the diffractive 
7r°-bremsstrahlung mechanisms as well as photon-photon fusion and photon-omega (omega- 
photon) exchange processes not discussed so far in the literature. The higher the energy, the 
two 7r°-bremsstrahlung contributions become better separated. The dotted line corresponds 



9 



to the photon-photon fusion mechanism. At the LHC energy and in the rapidity region 
—2 < y n o < 2 it even dominates over the diffractive mechanism. The cross section for the 
7r°-bremsstrahlung contribution at the LHC energy and at midrapidity is much smaller than 



e.g. for the production of heavy quarkonia: J /if) |42j, 
experimental measurement there would be a challenge. 



T |43| or Xco [35|, |44{. Clearly an 




FIG. 5. The distribution of vr in rapidity at ^/i = 45 GeV (ISR), 500 GeV (RHIC), and 14 TeV 
(LHC). The 7r°-bremsstrahlung contribution (black solid lines) and W7 (juj) exchanges (violet 
dashed lines) peaks at forward (backward) region of y n o , respectively, while 77 fusion (blue dotted 
lines) contributes at midrapidity. In this calculation we have used A^r = = 1 GeV of the 
hadronic form factors. No absorption effects are included here. 



Let us look now how absorption effects discussed in the theory section (see Figj2]) can 
modify the results obtained with the bare amplitudes (see FigfTJ). In FigJ2]we present, in 
addition, individual contributions for the 7r°-bremsstrahlung mechanism. We observe a large 
cancellation between the two terms in the amplitude [between the initial (p-exchange) and 
final state radiation (direct production)]. Because of destructive interference of bare and 
absorptive correction amplitudes, the resulting cross section is by a factor 2 to 3 smaller 
than that for the bare amplitude. The difference between the solid (An = = 1 GeV) and 
dashed (Ajv = 0.6 GeV and A w = 1 GeV) curves represents the uncertainties on the form 
factors. 

At large y^o another mechanism may come into the game - diffractive excitation of nucleon 
resonances. The resonances may occur when the energy in the irN subsystem W^n £ 7£> 
where 71 is the nucleon resonance domain. If (Wis) (y n o) G 1Z or (H23) iy-n ) £ 7Z then 
an extra strength due to resonance excitation may occur. In FigJT] we present the average 
value of subsystem energies (W13) and (W23) as a function of y T o at y/s = 0.5, 14 TeV. 
Only some nucleon resonances can be excited diffract ively. At the LHC they can occur for 
8 < \y n o\ < 11 and at the RHIC for 4.5 < \y n o\ < 7.5. One way to introduce resonances in 
the DHD model is to include them as intermediate states in the direct production term in 
Eq. (12.31) [see also FigJTfc)]. Calculating the contribution of diffractively produced resonances 
is more complicated and clearly goes beyond the scope of the present paper. The reader can 



find some theoretical attempts in Ref.|49 
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FIG. 6. The distribution of tt° in rapidity for y/s = 0.5, 14 TeV. In panel (a) we show individual 
contributions to the Born cross section. A large cancellation between the initial (p-exchange) and 
final state radiation (direct production) can be observed. In panels (b) and (c) the upper solid 
line corresponds to calculations without absorption effects, the lower solid line with absorption 
effects. The solid lines are for Ajv = = 1 GeV while the dashed lines are for A at = 0.6 GeV and 
A^ = 1 GeV. 



In FigJB]we show corresponding distribution in proton pseudorapidity again without (up- 
per lines) and with (lower lines) absorption effects and for two sets of A„- and An parameters. 
At the LHC protons could be measured by the ALFA (ATLAS) or TOTEM (CMS) detectors. 

The effect of absorption on transverse momentum spectra of protons and neutral pions 
is more complicated. In Figj9] we show distribution in transverse momentum of outgoing 
protons. Absorption causes a transverse momentum dependent damping of the cross section 
at small p± }P and an enhancement at large (compare upper and lower solid line). 

In Fig[l0]we show distribution in transverse momentum of 7T° meson. As in the previous 
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FIG. 8. The distribution in pseudorapidity of protons in the backward (left panel) and the forward 
(right panel) hemisphere at i/i = 14 TeV and for y n o > 0. Here = = 1 GeV (solid line) or 
A N = 0.6 GeV and A^ = 1 GeV (dashed line). 



figure we show results without and with absorption effects. The distributions are peaked at 
p ±t7T ~ 0.2 GeV. 

In Fig{TTJwe show distribution in the square of four-momentum transfer between initial 
and final protons. In panels (a) and (b) we show the separate contributions of different 
exchange terms. As in the previous figure we show results without and with absorption 
effects. One can observe much large tails of distributions in ti than in t 2 (y^o > was 
assumed) . 

In Figjl2]we show distribution in two-dimensional space £2) for the 7r°-bremsstrahlung 
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FIG. 9. The distribution of outgoing protons in transverse momentum for y/s = 14 TeV and for 
y n o > 0. As in the previous figure we show results without and with absorption effects. Here 
A at = = 1 GeV (solid line) or A^ = 0.6 GeV and A^ = 1 GeV (dashed line). 

contribution at y/s = 14 TeV (top panels) and y/s = 500 GeV (bottom panels) without (left 
panel) and with (right panel) absorption effects. The distributions in t\ or £2 are different 
because we have limited to the case of y n o > only. The distributions discussed here could 
in principle be obtained with the TOTEM detector at CMS to supplement the ZDC detector 
for the measurement of neutral pions. Similar analysis could be done by the ALFA detector 
for proton tagging at ATLAS. 

The pion energy spectrum for y n o > drops relatively slowly with pion energy which is 
shown in Fig{13l We show results without and with absorption effects. 

In FigJT3]we compare distribution in invariant mass of the forward produced pn° system 
for the 7r°-bremsstrahlung contribution and y n o > 0. The discussed here the pp — > ppn° 
process gives a sizable contribution to the low mass (Mx > fn P + "V ) single diffractive 
cross section. 

In Fig|T5] we show correlation function in azimuthal angle between outgoing protons. 
As can be seen in panel (a) the 7r°-bremsstrahlung contribution is peaked at back-to-back 
configuration (<pi2 = it). For comparison, the contribution for other mechanism 77 fusion and 
70; fusion are peaked at (fru = vr/2 and are much smaller. We observe a strong cancellation 
between the initial and the final state radiation. There [see panels (b) and (c)] is a sizable 
difference in shape between the result obtained with the bare amplitude and the result with 
inclusion of absorption effects. We doubt if such a correlation can be measured at the LHC 
in the future. 

In FigfTBlwe show distribution in two-dimensional space (t 2 ,M 13 ). One can observe dif- 
ferent behavior of slope in four-momentum transfer squared ti for different masses of the p7r° 
system. A similar effect was observed for pp — > p(n7t + ) [l8[ and np — > (pir^)p [2l[ reactions 
at much lower energies. As can be seen in FigsfIT] and O the large contribution comes from 
the 7r-exchange diagram and the baryon-exchange terms are suppressed due to amplitude 
cancellations. The differential cross section peaks for invariant masses close to threshold 
and disappears rapidly with increasing invariant mass, giving an approximately exponen- 
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FIG. 10. The distribution of ir° mesons in transverse momentum for y/s = 0.5, 14 TeV and for 
y n o > 0. In panel (a) we show individual contributions to the Born cross section. In the other 
panels as in the previous figures we show theoretical uncertainties (in e.g. form factors). Here 
A N = A n = 1 GeV (solid line) or = 0.6 GeV and A n = 1 GeV (dashed line). 



tial behavior for large masses. The absorptive effects could be partially responsible for the 
irregular structure in two-dimensional space fe, M i3 ) at small and M13 ~ 1.3 GeV. 

In FigJT7]we show corresponding two-dimensional distributions in (y n o,p tjn o) for y n o > 0. 
Sizable correlations between pion rapidity and transverse momentum can be observed which 
is partially due to interference of different components (amplitudes). 

In Table H] we have collected numerical values of the integrated cross section cr^^ after 
taking the forward region (y^o > 0) into account for exclusive production of ir° at different 
cm. energies \/s. Our results depend on the A parameters of the hadronic form factors. The 
cross section obtained from ISR experiments (see e.g., Ref . ["Lsj ) are roughly reproduced. 

For completeness, in FigJT8] (left panel) we compare the photon-odderon and odderon- 
photon contributions with the 77 contribution. We show results for B = 6 GeV -2 and two 
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FIG. 11. Distribution in the four- momentum transfer squared between initial and final protons at 
yfs = 14 TeV and for y n o > 0. In panels (a) and (b) we show individual contributions to the Born 
cross section. The theoretical uncertainties are shown in panels (c) and (d). Here A^r = A T = 1 GeV 
(solid line) or An = 0.6 GeV and A,,- = 1 GeV (dashed line). 



different estimates of the 7p — > 7t°p cross section (energy independent) as specified in the 
figure caption. The total cross section for the odderon contributions, corresponding to the 
HERA upper limit, is less than 20 nb in the rapidity region \y n o\ < 2.5. The corresponding 
curve is more than an order of magnitude larger than the photon-photon contribution. 

Finally, in FigJTH] (right panel) we make similar comparison of contributions of the two 
mechanisms for transverse momentum distribution of neutral pions for different slope pa- 
rameters and \y n o\ < 2.5. The curve corresponding to the HERA upper limit is considerably 
larger than the photon-photon contribution starting from p± <7r o > 0.2 GeV. Even with the 
Ewerz and Nachtmann limit, one can observe deviations from the 77 curve at transverse 
momenta p± t7T ° > 0.3 GeV. The cut on meson p_i_ i7r o should enhance relative odderon contribu- 
tion. In principle, the ALICE collaboration could try to measure the transverse momentum 
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FIG. 12. Distribution in (ti,t^) for the 7r°-bremsstrahlung contribution at y/s = 14 TeV (top 
panels) and y/s = 500 GeV (bottom panels) and for y n o > without (left panels) and with (right 
panels) absorption effects. Here A at = = 1 GeV. 

distribution of exclusively produced neutral pions. 

IV. A COMMENT ON SINGLE DIFFRACTIVE CROSS SECTION AT LOW PRO- 
TON EXCITATIONS 

The measurement of an inelastic proton-proton cross section is one of the standard and 
obligatory measurements at each collision energy. At the LHC single diffraction (SD) and 
double diffraction (DD) processes constitute a large contribution to the inelastic cross sec- 
tion (about a half). Unfortunately it is very difficult to truly measure the cross section 
for the low mass excitation at the LHC and often educated extrapolations are required. 
Usually 1/M 2 triple- Regge fit is used for this purpose. Do we have the expertise on the 
very low mass excitations? This issue was critically discussed recently j49|. The authors 
presented predictions of a dual-Regge model with a nonlinear proton Regge trajectory [|30[ 
with parameters fitted to the single diffractive cross section measured at low energies (for a 
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FIG. 13. Energy spectrum of pions at y/s = 0.5, 14 TeV and for y n o > 0. In panel (a) we show 
individual contributions to the Born cross section. Theoretical uncertainties are shown in panels 
(b) and (c). Here A N = A n = 1 GeV (solid line) or Aat = 0.6 GeV and = 1 GeV (dashed line). 




review of the low energy SD data see e.g. 22|, |5l|]). In their fit the low mass excitation is 



dominated by the excitation of the proton resonances iV*(1440) with J p = \ + and iV*(1680) 

with J p = | + . While the presence of the latter is rather natural - it is a member of the 
same Regge trajectory as proton - the huge contribution of the Roper resonance is not too 
clear to us. The low-energy experimental SD data [5l[ show up a huge peak at the nominal 
position of the Roper resonance. This is the region where the absorbed Drell-Hiida-Deck 
mechanism (the nonresonant background model) predicts an enhancement (see Fig JT4l . The 
arguments a gain st large Roper contribution in single diffraction at high energies were ex- 
posed in Ref . |52| . We wish to emphasize that the DHD contribution was not included in the 
analysis of the SD mass spectrum in j49j where only a purely mathematical fit was used. 
The fitted background seems to have quite different properties than the discussed here DHD 
mechanism with absorption (different both in Mx and in t). In our opinion inclusion of 
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FIG. 14. Distribution in proton-pion invariant mass M13 at ^/s = 0.5, 14 TeV and for y^o > 0. 
In panel (a) we show individual contributions to the Born cross section. Theoretical uncertainties 
are presented in panels (b) and (c). Here A at = A n = 1 GeV (solid line) or An = 0.6 GeV and 
A n = 1 GeV (dashed line). 



a realistic absorbed DHD contribution could dramatically change, or even eliminate, the 
contribution of the Roper resonance. This issue requires further studies. 



The resonances contributing to the SD cross section discussed in [49| naturally contribute 



also to the pp — > ppir channel and the corresponding cross section is 

= ^sd x BR(N* Vvr) x i . (4.1) 

The last factor comes from the fact that the considered diffractively excited baryon reso- 
nances have isospin I — |. The branching fractions BR(N* — > Ntt) have been measured [53[ 
and are about 65% for both discussed states. The same situation occurs in the pp — » p(n7r + ) 
and np — > (pn~)p reactions (a factor 2 larger cross section), where no clear signal of the 
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FIG. 15. The distribution in azimuthal angle between outgoing protons for y/s = 0.5, 14 TeV. In 
panel (a) we show individual contributions to the Born cross section. Theoretical uncertainties 
are presented in panels (b) and (c). Here Ajv = = 1 GeV (solid line) or A^r = 0.6 GeV and 
A n = 1 GeV (dashed line). 



Roper Af*(1440) resonance was identified (see e.g., 18N20|) while the A^*(1680) resonance 
was observed. □ The situation should be better clarified in the future. The discussed there 
resonances were not included in our analysis but could be included in principle. 
Our DHD mechanism contributes to the single diffraction cross section as 



a 



DHD 



'sd — 3^^p p7r o. (4.2) 

The factor 3 comes from the isospin symmetry of the NNn coupling constant. Taking our 
numbers from Table [J we predict certainly not a negligible contribution to the total inelastic 



2 In Ref.jlS] results on diffractive dissociation of protons into (nir + ) in pp collisions at the CERN ISR 
^/s = 45 GeV energy and the cross sections <J P p^.p(mr+) — (400 ± 110) /ib, T pp ^ pA r* (1680) — (170 ± 60) /xb 
have been reported. 
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(c) (d) 

FIG. 16. Distribution in fa, M13) for the 7r°-bremsstrahlung contribution at yfs = 14 TeV and for 
y^o > without [panel (a)] and with absorption effects in the final state only [panels (b) and (c)] 
and absorption effects in the initial state only [panel (d)]. Here Am = A n = 1 GeV [panels (a) and 
(b)] and Aat = 0.6 GeV, = 1 GeV [panels (c) and (d)]. 

cross section at high energies (for both-side SD the crf^^o should be multiplied by a factor 
2). To our knowledge the DHD contribution is not included in the existing Monte Carlo 
codes simulating high-energy diffractive processes. 

V. CONCLUSIONS AND DISCUSSION 

In the present analysis we have calculated differential cross sections for the exclusive 
pp — > pprc reaction at high energies relevant for RHIC and LHC We have included the 7r°- 
bremsstrahlung from the initial and final state, diffractive 7r°-rescattering, photon-photon 
fusion and photon-omega (omega-photon) fusion processes. The diffractive irN and NN 
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FIG. 17. Distribution in (y n o,p t7T o) at yfs = 14 TeV for the 7r°-bremsstrahlung contribution and 
y w a > without (left panel) and with (right panel) absorption effects. Here A^r = A^ = 1 GeV. 

TABLE I. The integrated value of cross sections in \ib for the pp — > ppifi reaction at yfs = 45 GeV 
(ISR), 500 GeV (RHIC), and 14 TeV (LHC) and y n o > is taken into account only. The lower 
limit corresponds to the result when Xjy = 0.6 GeV and X w = 1 GeV are imposed in calculations 
and the upper limit when A at = A^ = 1 GeV. 



Model 


y/s = 45 GeV 


^fs = 500 GeV 


y/s = 14 TeV 


No absorption 
Absorption in initial state 
Absorption in final state 


103 - 146 
46- 76 
60- 91 


177 - 251 
62 - 125 
84 - 139 


402 - 575 
94 - 357 
128 - 290 



rescattering amplitudes have been related to the total irN and NN cross sections. The 
Donnachie-Landshoff parametrization has been used for energy dependence of the latter. 
Absorptive effects have been included in addition. They lower the cross section by a factor 
2 to 3; see Tabled! 

We have found very large cross sections of the order of mb. The total (integrated over 
phase space) cross section is almost energy independent. The dominant contributions are 
placed at large rapidities. The larger cm. energy, the larger rapidities are populated. On 
the other hand, the diffractive contribution is absent at midrapidity (y n = 0). The higher 
the collision energy the larger the unpopulated region. This opens a window for other 
mechanisms with much smaller cross section. For example at the LHC the two-photon 
fusion mechanism "wins" with the diffractive mechanisms at y n ~ 0, where the diffractive 
contributions are very small. However, the transverse momenta of neutral pions in this 
region are very small and therefore such pions are very difficult to measure. The 70; or 
wy exchanges have been found to be significant only in backward or forward rapidities, 
respectively, and are small at midrapidities due to u reggezation. In principle, also a 2 - 
pomeron and pomeron-a 2 exchanges or p°-odderon and odderon-p° exchanges could play 
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FIG. 18. Rapidity distribution of neutral pions (left panel) produced in 70-fusion and 07-fusion 
(black upper solid lines) compared to the 77 contribution (blue lowest solid line) for ^fs = 14 TeV. 
Individual contributions of photon-odder on (short dashed line) and odderon-photon (long dashed 
line) are shown separately. We show predictions for the HERA upper limit (^■y P ^ 7T o p = 49 nb) 
and for the Ewerz-Nachtmann estimate {(T^ p ^^o p = 6 nb). In the right panel we make similar 
comparison of contributions of the two mechanisms for transverse momentum distribution of 7r°'s 
in the rapidity region —2.5 < y w o < 2.5. For the odderon contributions we have used different 
values of slope parameters B = 4, 6, 8 GeV~ 2 . 



some role but not at midrapidities. In addition, it is rather difficult to make for them 
realistic predictions. A larger cross section than predicted here at midrapidities would be 
an interesting surprise. 

We have shown several other differential distributions. If one limits to separate regions 
of y w o < or y^o > (one-side excitation), then the distributions in proton transverse 
momenta p\ t and p2t ar e quite different - one reflecting the pion/nucleon exchange and the 
second reflecting the pomeron exchange. The same is true for the t\ and t% (transferred 
four-momentum squares) distributions. Analysis of such details would be a useful test of 
the model. The distribution in the mass of the excited ir°p system peaks at small M np 
and quickly drops when the mass increases. Such a distribution reminds the spectral shape 
of the Roper resonance fitted recently to an old single- diffract ive data. We have obtained 
an interesting correlation between the mass of the excited system and the slope of the t 
distributions well represented in a two-dimensional plot d f° M (t, M). Similar effects were 
observed in the past for the pp — > p(nir + ) and np — > (p7r~)p reactions at the CERN ISR and 
Fermi National Accelerator Laboratory (Fermilab). 

Particularly interesting is the distribution in azimuthal angle between outgoing protons, 
not studied so far, including low-energy pp — > p(nn + ) and np — > {pit~)p reactions measured 
in the 1970's at ISR and Fermilab. The distribution has a maximum at relative angle 0x2 = ^ ■ 
The detailed shape of the distribution is, however, very sensitive to the relative contribution 
of different ingredients of the model. The sensitive nature of the cancellation between 
proton-exchange and direct production amplitudes leads to a situation where minor changes 
in the parametrizations of these amplitudes can have large effects on discussed distributions. 
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Experimental analysis of such a distribution would therefore help in fixing model parameters 
such as cutoff parameters of hadronic form factors, not known very precisely. 

At the LHC the tt° mesons could be measured with the help of zero degree calorimeters. 
Such measurements are possible only at rather large pseudorapidities \rj n o\ > 8 — 9. The 
pions at midrapidities have rather small transverse momenta so their registration is probably 
very difficult. On the other hand, protons could be measured with the ALFA detector of 
ATLAS or the TOTEM detector associated with the CMS main detector. The latter may 
be more difficult from an organization point of view. 

The reaction discussed here is interesting also in a much broader context. First of all, it 
may constitute a sizable fraction of the pion inclusive cross section at very forward/backward 
(pseudo)rapidities. A comparison with nondiffractive Monte Carlo code would be therefore 
very valuable. Second, it leads to a production of very energetic photons (~ 0.5 — 2 TeV) 
from the decay of the forward 7r°'s. These two issues will be a subject of future investigations. 
Third, the DHD mechanism sizably contributes to the single diffractive cross section and 
as a consequence to the total inelastic cross section. This contribution is not included in 
any existing Monte Carlo code. Needless to say, these codes are used when extrapolating 
the measured high-mass SD cross section down to the ttN threshold, which obviously leads 
to an underestimation of the extracted (measured and extrapolated) cross section for single 
diffraction and/or inelastic processes. Finally, because the cross section for the discussed 
reaction is large, detailed studies could help to test model(s) of soft absorption, so important 
in the context of more fundamental searches such as, e.g., exclusive production of the Higgs 
boson in diffractive processes. 

In the present paper we have calculated only contributions with intermediate protons in 
the ground state to the pp — > ppn reaction. There are also resonance contributions, due to 
diffractive excitation of some nucleon resonances and their subsequent decays into the p + n 
(p + 7t°) channels. The dominant contributions are due to N* resonant states being members 
of the nucleon trajectory. The iV*(1680) 5/2 + state is the best candidate. Although a huge 
contribution of the Roper resonance iV*(1440) was suggested recently j49|, as discussed in 
our paper, their contribution may be to some extent an artifact of a fit which does not 
include the DHD mechanism discussed in our paper, neither in the pp — > ppir nor in the 
pp — > pmr + channel. 

In the present analysis we have considered single exclusive 7r° production. One could 
think about immediate extension of the present study to double diffractive, double DHD 
mechanism producing two exclusive 7r°'s, not considered so far in the literature and not 
included in any Monte Carlo code. Again we expect a rather large cross section for such an 
inelastic process. This would be a competitive mechanism to the central 7r°7r° production 
via pomeron-pomeron fusion considered by us for 7r + 7r~ production; see [H, 0]. 

We have presented first estimates of the photon-odderon and odderon-photon contribu- 
tions based on the upper limit of the ^p — > ir°p cross section obtained at the HERA as well 
as estimates based on a nonperturbative approach of Ewerz and Nachtmann which makes 
use of chiral symmetry and PCAC. Based on the HERA upper limit we conclude that the 
cross section for the contribution to the pp — > ppn° reaction is smaller than 20 nb in the 
rapidity region |7/„.o| < 2.5. 

Any deviation from the 77 — > ir° contribution to transverse momentum distribution of 
neutral pions at midrapidity would be a potential signal of photon-odderon (odderon-photon) 
contributions. One can expect potential deviations from the photon-photon contribution at 
P±,tt° ~ 0.5 GeV. This requires dedicated studies if the considered process could be measured 
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by, e.g., the ALICE Collaboration at the LHC. 
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